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ABSTRACT 


In  this  thesis  a  reference-based  successive  interference  cancellation  (RSIC)  scheme  is 
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quadrature  phase  shift  keying  (QPSK)  modulation  in  orthogonal  frequency  division 
multiplexing  (OFDM)  system  are  analyzed.  Simulations  show  RSIC  schemes  give 
considerable  increase  in  performance  compared  to  a  system  where  interference 
cancellation  is  not  utilized. 
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I  INTRODUCTION 

A.  BRIEF  OVERVIEW 

Successive  interference  cancellation  (SIC)  is  a  well-known  technique  to  extract 
the  target  signal  from  multiple  signals.  The  multiple  signals  are  unwanted  interferences 
caused  by  signal  collision  or  overlapping.  SIC  uses  successive  signal  subtraction  using 
strong  reference  signals.  Conventional  SIC  technique  is  normally  applied  by  one  receiver 
intercepting  various  signals  from  multiple  transmitters.  Signal  of  interest  could  be 
strongly  contaminated,  i.e.  signal  strength  might  be  weak  compared  to  interference, 
which  causes  the  signal  of  interest  not  to  be  extracted  or  demodulated. 

In  certain  applications,  multiple  receivers  are  available.  A  new  scheme  called 
referenced-based  SIC  (RSIC)  is  proposed.  RSIC  successively  subtracts  interferences  by 
using  a  relatively  clean  signal  from  an  initial  receiver  in  a  multiplatform  receiver  system. 
By  deploying  a  receiver  in  a  favored  area,  an  initial  target  signal  can  be  chosen  which  is 
free  of  interference  (except  receiver  noise),  which  is  referred  to  as  reference  signal.  This 
signal  can  be  subtracted  in  a  receiver  placed  in  another  area  where  this  signal  is 
considered  the  only  interference  (again  except  for  its  own  receiver  noise).  Thus  the 
second  receiver  cancels  this  interference.  The  procedure  is  applied  to  the  rest  of  the 
platform.  Therefore  the  final  target  signal  can  be  extracted,  i.e.  performance  of  the  system 
is  greatly  improved. 

Two  different  reference-based  SIC  techniques  are  analyzed  in  this  thesis:  Non- 
Demodulated  RSIC  and  Demodulated  RSIC. 

B  SCOPE  AND  ORGANIZATION 

Two  signal  models  (for  Non-Demodulated  RSIC  and  Demodulated  RSIC  signal 
extraction  schemes)  are  proposed,  coded  and  simulated  using  Matlab  in  the  following 
chapters.  Three  important  variables  are  introduced:  the  amplitude  gain  of  an  interference 
signal  (a),  its  estimate  (/>'),  and  their  difference  (y).  In  order  to  illustrate  the  effect  of  the 
performance  of  one  variable  (say  y  ),  we  fix  the  value  of  that  variable  while  varying  the 
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another  (say  a).  In  this  work,  we  will  implement  this  approach  to  only  one  receiver  in  a 
multiple  receiver  scenario  (specifically  to  the  second  receiver  referred  to  as  Receiver  2). 

In  Chapter  II,  we  discuss  the  concept  of  a  reference-based  SIC  (both  Non- 
Demodulated  and  Demodulated).  We  discuss  how  to  configure  (by  positioning)  the  multi¬ 
receiver  platform  such  that  a  reference  signal  can  be  received  by  the  initializing  receiver. 

In  Chapter  III,  we  examine  the  proposed  scheme  by  simulating  multiple 
quadrature  phase  shift  keying  (QPSK)  receivers  using  Matlab.  Both  Non-Demodulated 
RSIC  and  Demodulated  RSIC  are  simulated. 

In  Chapter  VI,  we  examine  both  RSIC  schemes  with  multiple  quadrature  phase 
shift  keying  (in  OFDM)  receivers. 

In  Chapter  V,  we  conclude  the  thesis  and  provide  recommendations  for  future 
research. 
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II  REFERENCE-BASED  SUCCESSIVE  INTERFERENCE 

CANCELLATION  (RSIC) 

Receivers  pick  up  unwanted  signals,  noise,  and  interferences  from  targeting  areas. 
The  signal/s  of  interest  (SOI)  is/are  buried  among  noise  and  interferences.  Using 
estimated  reference  signals,  the  successive  interference  cancellation  (SIC)  technique  can 
be  applied  to  eliminate  the  interferences  and  extract  a  target  signal. 

The  conventional  SIC  works  with  one  receiver  [1].  In  a  multiple  receiver 
platform,  a  new  approach  called  reference-based  successive  interference  cancellation 
(RSIC)  is  proposed.  This  concept  involves:  deployment  of  multiple  receivers  in  favorable 
targeting  areas  and  then  the  receivers  subtract  interference  using  the  reference  signals 
passed  from  prior  receivers  until  SOI  (or  SOIs)  is  decoded. 

Two  reference-based  successive  interference  cancellation  techniques  (Non- 
Demodulated  and  Demodulated  RSIC)  are  addressed  in  this  chapter.  Obtaining  an  initial 
reference  signal  is  crucial  and  is  the  very  first  step  to  the  success  of  this  reference-based 
successive  interference  cancellation  process. 

A.  OBTAINING  THE  REFERENCE  SIGNAL 

For  the  RSIC  procedure  to  be  successful,  proper  configuration  or  positioning  of 
the  receivers  is  necessary. 

The  first  reference  signal  in  a  multiple  receiver  platform  is  crucial  to  the 
performance  of  the  whole  system.  Thus,  it  must  not  be  contaminated  with  the  signals 
from  the  other  emitters.  This  non-interfered  reference  signal  could  be  obtained  by 
positioning  the  receiver  such  that  other  interference  signals  are  not  present.  We  assume 
that  the  receiver  is  not  noise-limited,  i.e.  that  the  receiver  has  a  good  signal-to-noise  ratio 
(SNR). 

Conventional  cellular  system  employs  sectoring  antennas  to  reduce  co-channel 
interference  in  order  to  improve  link  reliability.  Shown  in  Figure  1,  four  base  stations 
employ  120°  antennas,  which  create  four  overlapping  areas.  By  positioning  the  receiver  in 
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Sector  A,  only  the  radio  signal  from  base  station  1  is  received,  and  it  is  used  as  the  initial 
reference  signal.  Radio  signals  from  both  base  station  1  and  2  are  received  by  Receiver  2 
in  Sector  B.  The  rest  of  the  receivers  can  be  positioned  in  Sector  C  and  D  to  intercept 
multiple  signals. 


Figure  1 .  Multiple  signals  received  in  various  sectors  (overlapping  antenna  coverage). 

B.  NON-DEMODULATED  RSIC 

Consider  the  case  of  deploying  four  receivers  in  four  targeting  areas  to  intercept 
all  of  the  four  base  station  radio  signals.  In  Figure  1,  a  signal  receiver  is  placed  in  each 
sector  (Sector  A,  B,  C  and  D)  ensuring  that  a  diversity  of  reception  is  achieved. 

According  Figure  1,  the  signals  received  from  various  sectors  are  listed  in  Table  1 
if  the  signals  are  of  equal  energy  (or  power). 
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Table  1.  List  of  receiving  signals  in  each  sector 


Sector 

Receivers 

Intercepted  signals 

A 

Receiver  1 

5)  + 

B 

Receiver  2 

+  s.,  +  n1 

C 

Receiver  3 

^i  +  ^2  ^3  ^3 

D 

Receiver  4 

^i  +  ^2  +  +  ^4  +  yi 4 

The  symbol  sm  represents  a  signal  from  base  station  m,  and  nn  is  the  addictive 
white  Gaussian  noise  (AWGN)  at  receiver  n. 

Clearly,  signal  2  ( s2 )  can  be  extracted  from  Receiver  2  by  subtracting  .v, .  An 
estimate  of  ^  can  be  sent  from  Receiver  1  to  Receiver  2.  The  signal  sent  from  Receiver  1 

to  Receiver  2  for  subtraction  is  called  a  reference  signal.  Repeating  the  same  procedure, 
signals  3  and  4  can  also  be  extracted  from  Receivers  3  and  4  by  using  the  subsequent 
reference  signal  passed  from  a  prior  receiver.  The  beauty  of  this  algorithm  is  that  the 
receivers  do  not  need  knowledge  of  the  information  content.  A  receiver  simply  passes  the 
reference  signal  to  subsequent  receivers  for  interference  cancellation. 

Assuming  equal  energy  for  all  signals  is  impractical.  Thus,  the  wireless 
transmission  amplitude  loss  or  gain  needs  to  be  considered  in  the  model.  Three 
interference-related  variables  are  used  to  model  received  signals.  Alpha,  am  n ,  is  the 

amplitude  gain  of  unwanted  signals  Sm  from  base  station  (emitter)  BSm  leaking  into 
receiver  Rn .  Beta,  pm  n ,  is  the  estimate  of  am  n  (amplitude  gain  of  reference  signal  Sm  ) 
passed  from  prior  receiver  Rm  to  subsequent  receiver  Rn  .  We  can  designate  ym  n  to  be 
the  difference  between  a  and/?,  i.e. 

/nun  -  am,n  ~  Pm,n  (2-1) 

Thus  the  received  signal  for  Receiver  1  is 


y\~s\  + 
5 


(2.2) 


where  the  received  signal,  yn ,  for  an  arbitrary  receiver  Rn  is: 


y2=s2  +  al2sl+n2 


(2.2) 


y?>  —  ^3  3^2  ^1  3^1  ^3 


(2.3) 


Jb  —  ^4  +  ^3,4^3  +  a2AS2  +  a\AS\  +  ^4 


(2.4) 


The  amplitude  gain,  ax  is  assumed  fixed  in  this  thesis.  In  practice,  it  may  have  to 

be  estimated  since  it  may  be  random  in  nature.  For  example,  it  may  be  distributed  from 
0.8  to  1.2,  where  the  0.4  deviation  maybe  due  to  proximity  of  receptors  in  relation  to 
emitter. 


If  we  let  Sj  =  sx  +  «,  be  the  estimate  of  desired  signal  in  Receiver  1,  then  ,v,  is  the 
initial  reference  signal  and  is  passed  from  Receiver  1  to  Receiver  2  for  cancellation 
procedure.  Receiver  2  receives  the  reference  signal  ,  including  its  own  AWGN  noise.  It 
performs  the  reference-based  interference  cancellation  process  as  dictated  by 

S2=y2-p\JX+fl2 

,  w  .  v  (2-5) 

=  {S2  +  «1,2J1  +n2)-  A, 2  (S1  +  ni  +  n2  ) 


Assuming  h2  is  negligible  (due  to  possibly  shorter  distance  between  two  receivers 
compared  to  distance  from  emitters)  then 


S2  5,2+(cTi  2  Pi'2 )  ‘h  ^2  P\.2^\ 

(2.6) 

y 1,2  —  a\,2  ~  P\,2 

(2.7) 

Thus, 

S2=S2+ri,2Sl+n2-fil,2ni 

(2.8) 

At  Receiver  3,  reference  signal  s3  can  be  calculated  via 

y 3  —  £3  +  (%2  3^2  ^1  3*^1  ^3 

(2.9) 

^3  —  53  —  P2,2p2  ~  A, 3*1 

(2.10) 
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S3  —  (s3  +  <^2,3‘s2  ^1,3^1  ;*3  )  Az.3  (^  "*"  Ti^l  ^2  P\,in\  )  A, 3  (‘h  ~*~  ^1  )  (2-11) 

? 


where  we  assume  the  noises  due  to  the  receipt  of  reference  signals  are  negligible.  Thus, 


A  ‘y3+(Q;2,3  /^2,3  )^2  (^1,3  A, 3  Pi,  l7x,  2)^1 
~  Pl,3n2  Pl,A.2n\  ~  P\,2n\  +n3 

(2.12) 

y 2,3 —  (^2,3 — Pi, 3 ) 

(2.13) 

11 

iT 

1 

> 

(2.14) 

Then, 

A  ^3  Yl,3S2  (7l,3  P 2,2,7 \, 2)^1  Pl,2ni  Pl,2,P\, inx  P\, 3^1  +  ^3 
—  53  + /2,3s  2  (^1,3  _  Pi,  3/l,2  )  ‘h  ~  Pl,3n2  (A, 3  A,  2  ~  P\,i)n\  +n3 

(2.15) 

With  Receiver  4,  reference  signal  ,v4  can  be  extracted  via 

T4  —  54  +  ^3,4^3  +  ^2,4^2  +  +  ,?4 

(2.16) 

A  —  T4 — A3, 4^3  ~  Pia^i  ~  P\aK 

(2.17) 

Thus, 

sA  =  (s4  +  a3  4s3  +  «2>452  +  +  n4 ) 

_  A3, 4  [^3  +  72,3s  2  {,7 1,2  ~  Pl,2,Y\,2  )51  ~  Pl,2n2  Jr{yPl,iP\,2  ^1,3  )  ,71  W3 
-  Pi  A  (52  +  Yx,lSX  +n2~  Px,2nx  )  -  A, 4  U  +  ) 

(2.18) 

Then, 

^4  =  ^4  +/3,4^3  +(/2,4  _  A3, 4^2, 3)  ^2  —  ^3,4  (7"l,3  ~  A2, 3/1,2  )  +  Pia7\,2  J  ‘h 

A3, 4^3  +  (PiaPi,!  ~  PlA  )  n2  [A2,4Al,2  —  Al,4  —  A3, 4  ^PiaPx,!  ~  PlA  )]  n\  +  W4 

(2.19) 

Non-Demodulated  RSIC  technique  does  not  require  the  receivers  to  demodulate 
the  intercepted  signals  to  perfonn  interference  cancellation.  The  information  of  the  signal 
content  is  not  required.  If  demodulation  is  of  interest  then  it  can  be  done  after  the  fact. 
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C.  DEMODULATED  RSIC 

In  developing  the  algorithm  of  reference-based  successive  interference 
cancellation  technique,  an  interesting  question  is  raised.  What  if  the  receiver  demodulates 
its  SOI  and  sends  the  demodulated  SOI  as  reference?  Is  there  a  perfonnance 
improvement?  In  the  demodulated  scheme  (now  appropriately  called  D-RSIC), 
demodulation  is  performed  in  each  receiver,  which  allows  receivers  to  make  decisions 
(detections)  on  intercepted  signals  based  on  the  modulation  technique  the  target  emitters 
are  using. 

As  in  Non-Demodulated  RSIC,  reference  signals  come  from  prior  receivers.  The 
major  difference  is  that  the  reference  signal  will  be  demodulated  and  re-modulated  before 
transmission  to  the  subsequent  receiver.  The  advantage  of  this  extra  step  mitigates  noise 
prior  to  transmission.  It  prevents  noise  and  error  accumulation  in  the  reference  signal. 
Subsequent  receivers,  especially  those  at  the  latter  part  procedure,  benefit  from  this  step. 

Adding  the  demodulation  procedure  into  Receiver  1,  the  detection  is  given  by: 

=  dec  (kj  +  nx )  (2.20) 

where  S1  is  the  new  initial  reference  signal  and  is  sent  from  Receiver  1  to  Receiver  2  and 
the  operation  “dec”  means  decide  which  symbol  is  received. 

At  Receiver  2,  received  reference  signal  1  is  used  for  interference  cancellation 
procedure.  Reference  signal  2,  s2  ,  can  be  calculated  by  demodulating  (2.8)  via 

s2  =  dec[y2 -/?12sj  +n2)  (2.21) 

The  same  demodulation  and  interference  cancellation  procedures  are  repeated  at 
subsequent  receivers.  Simulations  are  reported  in  Chapter  III,  where  significant 
improvement  in  receiver  performance  utilizing  this  technique  compared  to  Non- 
demodulated  RISC  is  demonstrated. 
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Ill  QPSK  SIMULATION 


A.  COMPARISON  OF  NON-RSIC,  NON-DEMODULATED,  AND 

DEMODULATED  RSIC  TECHNIQUE 

1.  Comparison  of  Three  Techniques  in  a  Two  Receiver  System 

Simulations  are  coded  in  Matlab  to  demonstrate  Non-RSIC,  Non-Demodulated 
RSIC,  and  Demodulated  RSIC  techniques  with  a  two-receiver  system  using  Matlab.  The 
signals  are  QPSK-modulated.  The  first  technique  (Non-RSIC)  is  one  without  benefit  of 
any  cancellation  schemes.  Receiver  2’s  received  signal  includes  signals  from  base 
stations  1,  2  and  noise.  The  Non-RSIC  scheme  demodulates  the  intercepted  signals 
without  the  benefit  of  RSIC.  In  the  second  technique  (ND-RSIC),  Receiver  2 
demodulates  the  intercepted  signals  by  subtracting  the  reference  signal  received  from 
Receiver  1.  In  third  technique  (D-RSIC),  reference  signal  1  is  demodulated  and  re¬ 
modulated  by  the  Receiver  1  prior  to  sending  to  Receiver  2.  With  this  extra  detection  step 
at  Receiver  1,  the  noise  of  intercepted  signal  1  is  mitigated  which  benefits  the  subsequent 
receiver  dramatically  (by  using  a  “cleaner”  reference  signal). 

In  Figure  2,  we  show  the  results  from  simulation,  which  show  bit  error  rate  (BER) 
performances  for  QPSK-modulated  signals  utilizing  the  three  techniques.  The  outer  curve, 
Non-RISC,  gives  us  a  baseline  for  comparison.  Without  any  interference  cancellation 
technique,  the  interference  and  noise  cause  substantial  performance  degradation.  The 
second  curve  is  the  BER  curve  for  ND-RISC  which  shows  substantial  improvement.  The 
benefit  of  D-RSIC  further  increases  the  performance.  This  is  the  curve  closest  to  the 
theoretical  QPSK  BER  curve. 

In  this  simulation,  the  interference  amplitude  gain  factor  al2  is  set  to  0.5  and  yl  2 
(ccj  2  minus  /?,  2 )  is  set  to  -0.05. 
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Figure  2.  Comparison  of  three  techniques  at  Receiver  2  with  two  receivers  where 

signals  are  QPSK-modulated 

2.  Comparison  of  Receiver  Performance  in  a  Four-Receiver  System 

a.  Noise  Accumulation  in  Subsequent  Receivers 

In  this  simulation,  four  receivers  are  considered.  With  ND-RSIC,  Receiver 
2  sends  non-demodulated  reference  signal  2  to  Receiver  3  for  SIC  procedure.  Receiver  3 
follows  the  same  procedure  where  the  non-demodulated  reference  signal  3  is  passed  to 
Receiver  4.  All  a  and  y  values  used  in  this  simulation  model  are  shown  in  Table  2. 
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Table  2.  Interference  amplitude  gain  factors  used  in  the  simulation 


Amplitude  gain  factors 

Value 

Note 

A, 2 

0.50 

al3 

0.90 

A, 4 

0.90 

A, 3 

0.80 

A, 4 

0.80 

A, 4 

0.70 

A, 2 

0.55 

A, 2  =  A, 2  ~Y\, 2 

A, 3 

0.95 

A, 4 

0.95 

A, 3 

0.90 

A, 4 

0.90 

A, 4 

0.90 

y  1,2 

-0.05 

Various  target  signals  broadcasted  from  multiple  base  stations  can  be  intercepted 
in  this  multiplatfonn  receiver  system.  A  signal  of  interest  (SOI)  can  be  extracted  from  a 
specific  receiver  according  to  Chapter  II.  It  may  be  that  the  perfonnance  is  worse  for 
latter  receivers  than  earlier  receivers  in  the  procedure.  But  this  really  depends  on  the 
a  s’,/?  s’  as  will  be  shown  later.  For  example,  if  the  broadcasted  signal  from  base  station 
4  is  the  SOI,  the  demodulated  signal  quality  (bit  error  rate,  BER)  at  Receiver  4  may  be 
the  worst  due  to  noise  accumulation  in  the  reference  signal  as  it  is  passed  along  from 
receiver  to  receiver. 

In  Figure  3,  we  show  the  results  of  four  Non-Demodulated  RSIC  (ND-RSIC) 
QPSK  receivers  (with  fixed  a  and  y  value).  The  ND-RSIC  BER  performance  of 
Receiver  2  is  better  than  that  of  a  Non-RSIC  Receiver  2.  This  is  due  to  the  cleaner 
reference  signal  from  Receiver  1  (but  without  demodulation).  The  performances  of  the 
subsequent  receivers  are  degraded  due  to  noise  accumulation  in  the  reference  signals.  The 
BER  performance  for  Receiver  3  and  Receiver  4  are  also  shown. 
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Figure  3.  Four  ND-RSIC  QPSK  receiver  system  (noise  accumulation  causes  signal 

degradation) 


b.  Adjusting  a  and  /?  Value  Enables  Subsequent  Receivers  to 
Outperform  Prior  ones 

We  understand  that  noise  accumulated  in  the  prior  reference  signals  causes 
performance  degradation  in  subsequent  receivers.  Flowever,  it  is  possible  for  a  latter 
receiver  to  outperform  a  prior  receiver.  The  performance  is  a  function  of  y  and  fi .  For 
example,  by  changing  fj,  4  value  from  0.90  to  0.65  (rest  of  the  variables  in  Table  2  are 
remaining  the  same),  Receiver  4  outperforms  Receiver  3.  The  result  is  shown  in  Figure  4. 

Thus,  the  performance  of  receivers  can  be  affected  by  the  estimates  of  the 
amplitude  gains.  In  general,  the  better  estimates  yield  performance  gain  (modest  or 
otherwise). 
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Figure  4.  Receiver  4  outperforms  Receiver  3  by  reducing  [3  value  from  0.95  to  0.65 


B.  COMPARISON  OF  NON-RSIC,  NON-DEMODULATED,  AND 

DEMODULATED  TECHNIQUE  AS  A  FUNCTION  OF  a 

To  further  understand  the  influence  of  interference  factor  a  andy  ,  various  a 
and  y  values  are  tested  in  the  simulation.  For  simplicity,  a  or  y  is  only  allowed  to  vary 
in  Receiver  2  of  the  multiplatform  receiver  system.  This  helps  to  observe  the  variation  of 
receiver  performance. 

In  this  section,  two-receiver  and  four-receiver  systems  with  various  a  scenarios 
are  coded  in  Matlab  for  simulation. 
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1. 


Two  Receivers 


A  two  QPSK-receiver  system  is  tested  with  the  three  techniques  as  in  previous 
experiments,  but  with  various  a  .  Thus,  al2  is  the  only  variable  that  varies  in  this 


simulation. 


Three  RSIC  in  two  QPSK  receivers  system 


Three  RSIC  in  two  QPSK  receivers  system 


Figure  5.  Comparison  of  three  techniques  in  a  two  QPSK  receiver  system  with  various 

a  and  fixed  y 


Shown  in  Figure  5,  all  three  techniques  are  tested  with  four  ain  values  which  is 

increased  from  0.5  to  2  with  0.5  increments.  The  outer  curve  (labeled  Non-RSIC)  shows 
that  without  the  benefit  of  RISC,  the  BER  performances  suffer  greatly  when  ,  is 

increased.  For  example,  with  Non-RISD  it  would  be  difficult  to  extract  SOI  when  ,  is 
greater  than  1 .  With  ND-RSIC  technique,  the  performance  of  the  receiver  starts  to  suffer 
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as  «, ,  values  increase  despite  corresponding  increase  in  the  estimates  (recall  that  the  y  is 

fixed).  However,  D-RSIC  receivers  perforin  extremely  well  regardless  the  increasing  of 
a,  2  value.  In  conclusion,  BER  improves  dramatically  with  the  benefits  of  RSIC 

techniques. 


2.  Four  Receivers 


In  Figure  6,  we  show  the  performance  of  four  ND-RSIC  receivers  demodulating 
target  signals  of  interest.  The  four  plots  are  produced  using  the  values  of  variables  in 
Table  2,  except  ax ,  at  Receiver  2  is  changed  from  0.5  to  2.0  with  0.5  increments.  Clearly 

it  can  be  seen  that  the  performance  of  Receiver  2  is  greatly  influenced  by  increasing  ax , 

value  despite  corresponding  increase  in  the  estimates.  The  BER  performances  suffer 
greatly  as  a  function  of  increasing  a,,. 


Four  ND-RSIC  QPSK  receivers  demodulating  signals 


Figure  6.  BER  performances  of  a  four  ND-RSIC  QPSK  receiver  system  with  various  a 
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Despite  the  impact  of  the  increasing  a ,  2  at  Receiver  2,  the  performances  of  the 

subsequent  receivers  (Receivers  3  and  4)  are  minimally  affected  since  their  amplitude 
gains  and  corresponding  estimates  are  not  allowed  to  vary  in  the  simulation. 

To  combat  the  impact  of  the  increasing  alpha  on  receiver  performance,  the 
Demodulated  RSIC  technique  is  implemented  in  Matlab  involving  four  receivers.  The 
simulation  results  are  shown  in  Figure  7.  Notice  the  dramatic  improvement  regardless  of 
the  al2  value. 


Four  D-RSIC  QPSK  receivers  demodulating  signals 


Figure  7.  BER  performances  of  four  D-RSIC  QPSK  receiver  system  with  various  a 

values 


The  D-RSIC  BER  improvements  for  Receiver  3  and  Receiver  4  (compared  to 
ND-RSIC  Receiver  3  and  Receiver  4  BERs)  are  worth  mentioning.  First,  the 
performances  of  both  receivers  greatly  improve  (compared  to  ND-RSIC  performances) 
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since  accumulated  noises  in  the  reference  signals  are  mitigated.  Receiver  3’s  BER 
benefits  more  and  performs  as  well  as  Receiver  2.  Second,  Receiver  4  no  longer 
outperforms  Receiver  3  as  shown  in  the  ND-RSIC  simulation  (Figure  6).  The  influence  of 
the  amplitude  gains  is  mitigated  when  the  demodulation  procedure  is  added.  This  is 
another  advantage  of  the  D-RSIC  technique. 

In  retrospect,  there  could  be  a  concern  regarding  the  addition  of  the  demodulation 
procedure  to  the  receiver.  It  may  worsen  the  estimate  of  the  reference  signal  due  to 
decision  making  error.  Once  the  error  is  made,  the  original  transmitted  signal  is  not  fully 
recovered  (re-created)  due  to  few  eventual  errors  contributed  by  the  random  nature  of 
receiver  noise.  But,  the  previous  simulation  results  prove  the  benefit  of  actually  making 
decisions  to  mitigate  the  AWGN  noise. 

Another  benefit  of  adding  demodulation  procedure  is  shown  in  Figure  8.  The  plot 
on  the  left  shows  ND-RSIC  QPSK  Receiver  2  demodulating  signals  with  four  different 
a  1 2  values.  The  BER  plot  on  the  right  (using  the  same  a  values)  corresponds  to  where 

D-RSIC  technique  is  implemented.  The  effect  of  noise  accumulation  in  the  reference 
signals  is  nullified  by  D-RSIC  technique  regardless  of  increasing  ,  values. 


D-RSIC  QPSK  Receiver  2:  Demodulating  signals  with  various  alpha 


Figure  8.  Comparison  of  ND-RSIC  and  D-RSIC  QPSK  Receiver  2  with  various  a 

values 
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C.  COMPARISON  OF  NON-RSIC,  NON-DEMODULATED,  AND 
DEMODULATED  TECHNIQUE  AS  A  FUNCTION  OF  y 

In  this  section,  the  influence  of  interference  factor  y  is  analyzed.  In  the 
simulation  the  y  values  are  varied  to  investigate  the  effects  on  BER. 


1.  Two  Receivers 


In  this  simulation,  y|  2  is  the  variable  allowed  to  vary.  Since  yx ,  is  the  difference 
between  ax2  and  j3x2  (yX2  =  aX2  -  J3l2),  changing  yl2  value  (with  fixed  «12  scenario)  is 
literally  changing  [j]  2  value  as  well. 


Figure  9.  Comparison  of  three  techniques  in  a  two  QPSK  receiver  system  with  fixed  a 

and  various  y  values 


In  Figure  9,  the  comparison  of  three  techniques  with  various  yX2  in  a  two  QPSK- 
receiver  system  is  tested.  The  value  of  yxl  is  varied  from  -0.1  to  -0.025  with  0.025 
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increments.  The  deviations  between  these  four  BER  plots  are  very  small.  The 
performances  of  two  RSIC  techniques  are  slightly  improved  when  the  absolute  value  of 
1/  2  is  decreased.  The  factor  y  2  is  the  difference  between  amplitude  gain  and  its  estimate 

of  reference  signal  1  as  received  by  Receiver  2.  With  fixed  «12  value,  /?12  value  is 
increased  and  the  absolute  value  of  yl  2  is  decreased.  The  result  is  modest  improvement  in 
BER  as  the  estimates  get  closer  to  amplitude  gains. 

The  influence  of  changing  y  on  ND-RSIC  receivers  seems  to  be  not  as  significant 
when  compared  to  changing  a  value  (with  fixed  y  ).  Note  however  that  the  y  variation 
is  smaller  compare  to  a  variation  used  in  previous  experiment. 

2.  Four  Receivers 

The  simulation  results  of  a  four  ND-RSIC  QPSK-receiver  system  demodulating 
signals  with  various  yl2  are  shown  in  Figure  10.  The  yl2  value  is  varied  from  -0.1  to  - 

0.025  with  0.025  increments.  The  simulation  results  show  minor  variation  with  various 
yl2  values.  Receiver  3  is  outperforms  Receiver  4  due  to  specially  assigned  fj,  4  value 

according  to  Table  2. 

The  advantage  of  D-RSIC  technique  in  four-receiver  system  is  shown  in  Figure 
1 1 .  The  performance  of  the  receiver  BERs  is  increased  compared  to  ND-RSIC  BERs  in 
Figure  10.  The  performance  impact  of  improving  y12  value  is  modest.  Again,  the  ;/ 

variation  is  smaller  compared  to  a  variation  used  in  a  previous  experiment. 
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BER  BER 


Four  ND-RSIC  QPSK  Receivers  Demodulating  Signals 


Figure  10. 


BER  performances  of  a  four  ND-RSIC  QPSK  receiver  system  with  various  y 


Figure  1 1 . 


BER  performances  of  a  four  D-RSIC  QPSK  receiver  system  with  various  y 
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The  y  value  is  the  measurement  of  how  close  the  jB  approaches  to  a  value. 
Clearly,  closer  the  two  values  are  in  terms  of  absolute  value  is  the  better.  Theoretically, 
the  interference  can  be  totally  eliminated  if  y  =  0  (a  =  ft). 

Shown  in  Figure  12,  BER  performance  of  Receiver  2  is  improved  slightly  when 
yx  2  (which  is  negative  in  value)  increases.  Flowever,  it  decreases  in  absolute  value.  When 

yl2  value  increases  (positive  values  as  shown  in  Figure  13),  the  performance  worsens  for 
both  schemes.  Flowever,  D-RSIC  scheme  is  clearly  more  robust  than  ND-RSIC  scheme. 

These  two  scenarios  illustrate  that  the  lower  absolute  y  ,  value  results  in  the  BER 
performance  gain. 


Figure  12.  Comparison  of  ND-RSIC  and  D-RSIC  QPSK  Receiver  2  with  various  y 

(negative  y  values) 
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Figure  13.  Comparison  of  ND-RSIC  and  D-RSIC  QPSK  Receiver  2  with  various  y 

(positive  y  values) 


In  Figure  13,  the  BER  variation  is  very  small  when  yX2  value  is  below  0.1.  D- 
RSIC  receiver  has  better  immunity  to  changing  yx  2  values.  As  mentioned  before,  overall 
D-RSIC  BER  performances  are  better  than  ND-RSIC  BER  performances. 

The  simulations  in  this  chapter  prove  that  RSIC  techniques  boost  receiver 
performance  in  multi-platform  systems.  The  interference  factor  a  seems  to  significant 
degradation  on  receiver  performance  compared  to  that  of  the  y  factor.  But  we  note 
however  that  the  y  variation  is  smaller  compared  to  a  variation  used. 

D-RSIC  is  shown  to  be  the  better  cancellation  technique  compared  to  ND-RSIC. 
This  multi-platform  receiver  system  can  benefit  applications  designed  to  extract  multiple 
signals  from  target  areas. 


22 


IV  OFDM  SIMULATION 


A.  COMPARISON  OF  NON-RSIC,  NON-DEMODULATED,  AND 

DEMODULATED  TECHNIQUE  AS  A  FUNCTION  OF  a 

In  this  chapter,  the  RSIC  techniques  are  tested  with  the  orthogonal  frequency 
division  multiplexing  (OFDM)  technique  which  is  now  widely  used  in  many  wireless 
systems  including  the  popular  IEEE  802.1  la,g  WLAN  [2,3].  Simulations  utilizing  OFDM 
with  64  QPSK  modulation  subcarriers  with  6  symbols  cyclic  prefix  (CP)  are  coded  using 
Matlab. 

The  base  simulation  model  uses  the  same  scenario  discussed  in  the  previous 
chapter,  only  one  receiver’s  (Receiver  2)  variable  ( a  or  y  )  is  allowed  to  vary  at  a  time 
during  the  simulation.  The  values  for  these  variables  are  configured  using  Table  2. 

1.  Two  Receivers 

In  Figure  14,  we  show  Receiver  2  BER  results  with  the  three  previous  techniques 
in  a  two-OFDM  receiver  system.  The  ax  2  value  is  varied  from  0.5  to  2  with  0.5 

increments.  Despite  the  additional  complexity  of  the  OFDM  modulation,  the  BER  results 
are  almost  the  same  as  the  perfonnance  shown  in  Chapter  IV  using  QPSK  receivers.  This 
is  due  to  the  orthogonal  nature  of  the  modulation.  With  non-RSIC,  the  BER  suffers 
performance  degradation  with  increasing  ax  2  value.  With  ND-RSIC,  BER  performance  is 

better  compared  to  not  doing  any  cancellation  (Non-RSIC).  The  D-RSIC  technique 
provides  the  best  performance. 
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Figure  14.  Comparison  of  three  techniques  in  a  two-OFDM  receiver  system  with  various 

a  and  fixed  y  values 


The  BER  performance  of  OFDM  in  AWGN  is  the  same  as  the  BER  performance  of  the 
underlying  modulation  scheme  (QPSK  in  this  case)  due  to  orthogonality  of  the 
subcarriers  as  mentioned  before. 

Cyclic  prefix  (CP)  is  a  guard  time  added  to  the  start  of  an  OFDM  symbol  to 
eliminate  intersymbol  interference  (ISI)  [4].  The  drawback  of  adding  CP  is  a  loss  in 
throughput.  In  our  simulation,  the  number  of  samples  is  increased  from  64  to  80,  which 
also  increases  symbol  energy.  Thus  the  signal  level  needs  to  be  scaled  down  by  a  factor 
of  V64/80  to  keep  the  total  symbol  energy  at  the  same  level.  This  causes  nearly  a  ldB 
loss  in  BER  performance  as  shown  in  this  chapter. 
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2. 


Four  Receivers 


In  a  system  of  four  OFDM  receivers  the  value  of  a,  2  is  increased  from  0.5  to  2  in 

0.5  increments.  Both  ND-RSIC  and  D-RSIC  techniques  are  tested.  The  simulation  results 
are  shown  in  Figures  15,  16  and  17.  By  demodulating  the  reference  signals,  the  BER 
performances  of  D-RSIC  technique  are  better  than  that  of  ND-RISC  technique. 


Four  ND-RSIC  OFDM  Receivers  Demodulating  Signals  Four  ND-RSIC  OFDM  Receivers  Demodulating  Signals 


Figure  15.  BER  performances  of  a  four  ND-RSIC  OFDM  receiver  system  with  various 

a  values 
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Four  D-RSIC  OFDM  Receivers  Demodulating  Signals 


Figure  16.  BER  performances  of  a  four  ND-RSIC  OFDM  receiver  system  with  various 

a  values 


D-RSIC  OFDM  Receiver  2  Demodulating  Signals  With  Various  Alpha 


Comparison  of  ND-RSIC  and  D-RSIC  OFDM  Receiver  2  with  various  a 
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Figure  17. 


B.  DEMODULATED  TECHNIQUE  AS  A  FUNCTION  OF  y 


In  this  section,  the  influence  of  difference  factor  y  is  analyzed.  According  to  the 
QPSK  receiver  simulation  in  Chapter  IV,  the  influence  of  the  difference  factor  y  is  not 
as  dramatic  as  varying  amplitude  gain  factor  a  although  we  note  that  the  y  variation 
used  is  smaller  compared  to  a  variation  used  in  the  experiments  in  Chapter  IV. 


1.  Two  Receivers 

In  Figure  18,  we  show  the  comparison  of  three  techniques  in  a  two  OFDM 
receiver  system  with  fixed  «12  value  (0.5)  and  various  y]  2  values.  The  value  of  yl2  are 

negative.  It  is  increased  from  -0.1  to  -0.025  with  0.025  increments  (decrease  in  absolute 
values). 


Figure  18.  Comparison  of  three  techniques  in  a  two  OFDM  receiver  system  with  fixed  a 

and  various  y  values 
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BER  BER 


2. 


Four  Receivers 


In  Figures  19  and  20,  we  show  the  BER  results  of  two  cancellation  techniques.  As 
now  expected,  D-RSIC  technique  is  superior  to  ND-RISC  in  terms  of  performance. 


Figure  19.  BER  performances  of  a  four  ND-RSIC  OFDM  receiver  system  with  y 

various  value 
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Figure  20. 


BER  performances  of  a  four  D-RSIC  OFDM  receiver  system  with  various  y 

value 


In  Figures  21  and  22,  we  show  the  BER  curves  fir  Receiver  2  with  various  yl2  values. 
The  smaller  y  absolute  value  means  BER  performance  improves  correspondingly. 
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BER  C  BER 


2 1 .  Comparison  of  ND-RSIC  and  D-RSIC  OFDM  Receiver  2  with  various  y 

(negative  y  values) 


Figure  22.  Comparison  of  ND-RSIC  and  D-RSIC  OFDM  Receiver  2  with  various  y 

(positive  y  values) 
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In  this  chapter,  we  consider  OFDM  with  QPSK  modulation.  Both  RISC 
techniques  improve  the  BER  perfonnances  of  multi-platform  receivers.  The  D-RSIC 
technique  performs  the  best.  The  interference  factor  a  causes  significant  impact  on 
receiver.  With  D-RSIC,  the  perfonnance  is  not  as  sensitive  to  a  variations,  which  is  the 
great  benefit  of  this  technique. 
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V  CONCLUSIONS 


A.  SUMMARY 

Reference-based  successive  interference  cancellation  techniques  (RSIC) 
techniques  boost  the  BER  performances  of  multi-platform  receivers.  The  first  intercepted 
signal  (used  as  initial  reference  signal)  is  sent  to  the  subsequent  receiver  to  be  subtracted. 
Once  the  interference  is  eliminated,  a  new  reference  signal  is  sent  to  the  next  receiver. 
After  the  final  receiver,  the  final  target  signal  or  signal  of  interest  can  be  extracted. 

One  crucial  key  step  in  the  RSIC  procedure  is  to  obtain  an  interference-free 
reference  signal.  This  can  be  achieved  by  positioning  the  receiver  at  a  favorable  location. 
The  same  goes  for  other  receivers. 

Three  variables,  gain  factor,  its  estimate,  and  difference  gain  factor,  (a  J3 ,  and 
y  )  are  utilized  to  model  received  signals.  The  amplitude  gain  a  is  the  gain  of  unwanted 
signal  leaking  into  a  receiver.  The  variable  /?  is  the  estimate  of  a  for  that  receiver.  The 
variable  y  is  the  difference  between  a  and  ft . 

In  our  simulations,  the  perfonnance  of  a  receiver  is  greatly  influenced  by 
changing  a  .  Similarly  the  performance  is  also  influenced  by  varying  y  .  By 
implementing  the  D-RSIC  technique,  significant  improvements  in  BER  are  realized  due 
to  the  elimination  of  the  noise  accumulation  in  the  reference  signals.  The  D-RSIC 
technique  adds  demodulation  to  each  receiver  thereby  requiring  more  signal  processing, 
software  and  processing  time  compared  to  ND-RSIC  receivers,  but  the  benefits  of  the 
extra  step  pays  off  tremendously  in  BER  performance. 

B.  CONCLUSION 

The  performance  of  new  proposed  reference-based  successive  interference 
cancellation  techniques  are  validated  by  Matlab  simulations.  The  BERs  are  improved. 
Thus,  target  SOIs  can  be  extracted. 
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The  D-RSIC  technique  has  more  signal  processing  needs  than  ND-RSIC  but 
yields  a  tremendous  gain  in  BER  perfonnance. 

C.  FUTURE  WORK 

We  recommend  future  studies  to  explore  actual  estimation  of  a  with  f3 . 
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APPENDIX  A 


Figures  2  to  8  (QPSK  simulation  with  various  a  )  and  Figures  14  to  17  (OFDM 
simulation  with  various  a  )  are  generated  by  the  following  Matlab  codes: 

9-  •k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 
o 

%  Description: 

%  1.  Four  QPSK/OFDM  emitters  and  receivers 

%  2.  Three  techniques:  non-RSIC,  Non-Demodulated  and  Demodulated  RSIC 
%  3.  At  Receiver  2,  alphal2  changes  from  0.5  to  2,  Gammal2  =  - 
0.05 (fixed) 

%  4.  For  Thesis  Chapter  III  QPSK  and  Chapter  IV  OFDM  various  alpha 

9-~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k 
o 

clc;  clear  all;  close  all 
format  long; 

%  choose  simulation  model:  (l)QPSK;  (2) OFDM 

disp('RSIC  simulation  with  various  alpha  and  fixed  gamma  value') 
ch  =  input (' Choose  (l)QPSK  Simulation  or  (2) OFDM  Simulation  :  '); 

tic  %  start  clock  to  measure  performance 


%  Bit  length 


if  ch  ==  1 


QPSK  simulaiton 


bit  length  =  le6; 
elseif  ch  ==  2  %  OFDM  simulation 
%  Bit  Length 
M=4  ; 

k=log2 (M) ; 
n  fft=64; 
n  cyc=16; 


Number  of  symbols  in  constellation  (QPSK 
Bits  per  symbol 

Size  of  FFT (number  of  OFDM  carriers) 
Length  of  cyclic  prefix  (OFDM) 


=  4) 


n  symbols  =  le4; 


Number  of  symbols;n  symbols= (bit_count/k) /n  fft 


bit_length=k*n  fft*n  symbols;  %bl= (bits/symbol) * (FFTsize) * (symbols) 


end 

%  Range  of  SNR(dB)  to  simulate 
EbNo  =  -3:  1:  10; 

%  Convert  Eb/No  values  to  channel  SNR 
SNR  =  EbNo  +  10*1 og 10 (2)  ; 


%%  Define  variables  -  interference  amplitude  gain  factor 
%  Alpha:  Interference  signal  amplitude  received  at  receiver 


alphal3  = 

.9;  % 

Signal 

1 

leaks 

into 

Receiver 

3 

alpha23  = 

.9;  % 

Signal 

2 

leaks 

into 

Receiver 

3 

alphal4  = 

.8;  % 

Signal 

1 

leaks 

into 

Receiver 

4 

alpha24  = 

.8;  % 

Signal 

2 

leaks 

into 

Receiver 

4 

alpha34  = 

.7;  % 

Signal 

3 

leaks 

into 

Receiver 

4 

%  Beta:  Reference  signal  amplitude  received  at  subsequent  receiver 


betal3  = 

LO 

Oh 

%Receiver 

1 

to 

receiver 

3 

betal4  = 

.95; 

%Receiver 

1 

to 

receiver 

4 

beta23  = 

.9; 

%Receiver 

2 

to 

receiver 

3 

beta24  = 

.9; 

%Receiver 

2 

to 

receiver 

4 

beta34  = 

.  65; 

%Receiver 

3 

to 

receiver 

4 

%  Various  alphal2  and  fixed  gamma 

alphal2  =  0.5:  0.5  :2;  %  Change  from  0.5  to  2  with  0.5  increments 
gamma!2  =  -0.05;  %  Fixed  Gamma  value  (Gamma  =  Alpha  -  Beta) 
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%%  create  BER  buffers  (to  increase  Matlab  speed,  not  necessary) 

BER  =  zeros ( 4 , length ( SNR) ) ;  %  4:  FOUR  receivers 

%save  Rcvr2  SNR  of  various  alphal2 

BER  ND  RSIC  Rcv2  alpha=  zeros (length (alpha!2) , length (SNR) ) ; 

BER  D  RSIC  Rcv2  alpha  =  zeros (length (alphal2) , length (SNR) ) ; 

%%  Start  various  Alphal2  calculation  loop 

for  alpha  count  =1:  1:  length (alphal2 )  %  Calculate  each  alphal2 

betal2  =  alphal2 (alpha  count)  -  gammal2;  %  Changes  with  Alphal2 
for  RSIC  count  =1  :  1  : 3  %  Do  three  RSIC  techniques  one  at  a  time 
%1:  Non-RSIC;  2:  Non-demodulated  RSIC;  3:  Demodulated  RSIC 
%%  Start  the  main  calculation  loop  for  each  SNR 
for  snr  count  =1:  1:  length (SNR) 

%  Initiate  values  for  error  buffer 

Tx  Errors  =  [0;  0;  0;  0];  %  FOUR  rows  =  FOUR  receivers 
Tx_bits  =  [0;  0;  0;  0] ;  %  FOUR  rows  =  FOUR  receivers 
%  Keep  running  until  we  get  100  errors  in  Signal  1 
while  Tx  Errors  (1,:)  <  100 

%%  Generate  information  bits.  Four  rows  for  Four 
signals (emitters) 

uncoded  bits  =  round (rand ( 4 , bit  length));  %  4:  FOUR  receivers 
%%  Split  each  signal  into  two  streams,  for  Quadrature  Carriers 
I  =  uncoded_bits ( : , 1 : 2 : end) ;  %  get  odd  bits 

Q  =  uncoded_Joits  ( : ,  2  :  2  :  end)  ;  %  get  even  bits 

%  QPSK  modulator  set  to  pi/4  radians  constellation (Transmitter) 
Tx  sig  = . . . 

(  (I==0)  . * (Q==0) * (exp (li*pi/4) )  .  . . 

+  (1=0)  .  *  (Q==l )  *  (exp  (3*li*pi/4)  )  .  .  . 

+  (1=1)  .  *  (Q=l)  *  (exp  (5*li*pi/4)  )  .  .  . 

+  (1  =  1)  .  *  (Q=0)  *  (exp  (7*li*pi/4)  )  )  ; 
if  ch  ==  2  %  OFDM  simulation 
%%  OFDM  Modulation 
for  n  =  1  :  n  symbols; 

s  of dm  =  sqrt (n  f f t) *if f t (Tx  sig ( : , (n-1 ) *n  f ft+1 :  . . . 

n*n  fft) ' ,n  fft) ;  %Normalized  IFFT  operation 
s  ofdm  eye  =  sqrt (n  fft/ (n  fft+n  eye))  *  [s  ofdm(n  fft... 

-n  cyc+1  :  n  fft, :);  s  ofdm];  %Appending  cyclic  prefix 
OFDM  tx ( : , (n-1) * (n  fft+n  eye) +1  :  n* (n  fft+n  eye) ) . . . 

=  s  ofdm  eye';  %  Combine  subcarriers  into  one  to  tx 

end 

Tx_sig  =  OFDM_tx; 

end 


%%  Gaussian  Link  Noise 

No  =  1/10A (SNR (snr  count) /10);  %Noise  variance 

Noise  =  sqrt (No/2 )*.. .  %  4  ramdom  noises  for  four  channels 

(randn (4, length (Tx  sig) ) +li*randn (4, length (Tx  sig))); 
%%  Receiver  intercepted  signals 

%  including  target  SOI  and  interferences  from  adjeent  emitters 
Rx  received  =  zeros (4, length (Tx  sig));  %  Create  buffer 
%  Receiver  1:  Signal  1  +  Noise  1 
Rx_received ( 1 , ; )  =  Tx_sig(l,:)+  Noise(l,:); 

%  Receiver  2:  Signal  2  +  (alphal2  *  SI)  +  Noise  2 
Rx  received (2 ,: )  =  Tx  sig (2, :)+... 

(alphal2 (alpha_count)  *  Tx_sig ( 1 , : ) ) +  Noise(2,:); 

%Receiver  3;  Signal  3  +(alpha23  *  S2)  +(alphal3  *  SI)  +Noise  3 
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Rx_received ( 3 , : )  =  Tx_sig(3,:)+  (alpha23  *  Tx_sig (2, :))+.. . 

(alphal3  *  Tx_sig ( 1 , : ) ) +  Noise(3,:); 

%  Receiver  4:  S4  +  (a34  *  S3)  +  (a24  *  S2)  +(al4  *  SI)  +Noise  4 
Rx_received ( 4 , : )  =  Tx_sig(4,:)+  (alpha34  *  Tx_sig (3, :))+.. . 

(alpha24  *  Tx  sig (2, : ) ) +  (alphal4  *  Tx  sig (1, : ) ) +  Noise(4,:); 

0,0, _ 

o  o 

%%  Reference-based  successive  interference  cancellation 
procedures 

%  run  one  RSIC  technique  at  a  time 
%%  Non-RSIC 
if  RSIC_count  ==  1 

Rx  proced  =  Rx  received;  %  no  cancellation  procedure 

%%  Non-Demodulated  RSIC 
elseif  RSIC  count  ==  2 

%  Signal  1  is  used  as  reference  directly;  Intercepted  signals 
%  will  be  subtract  by  reference  signal  with  a  factor  Beta 
Rx  proced  =  Rx  received;  %  Create  buffer  and  save  Signal  1 
Rx_proced (2 ,  : )  =  Rx_received (2 ,  : ) - (betal2  *  Rx_proced ( 1 ,  : ) ) ; %S2 
Rx_proced ( 3 , : )  =  Rx_received ( 3 , : )  -  (beta23*Rx_proced (2 ,:)).. . 

-  (betal3  *  Rx_proced ( 1 , : ) ) ;  %  get  Signal  3  only 

Rx_proced ( 4 , : )  =  Rx_received ( 4 , : )  -  (beta34 *Rx_proced ( 3 ,:)).. . 

-  (beta24  *  Rx_proced (2 , : ) )  -  (betal4*Rx_proced ( 1 , : ) ) ;  %  S4 
%%  Demodulated  RSIC 

elseif  RSIC  count  ==  3 

%  Demodulate  1st  received  QPSK  signal (with  noise)  at  Receiver  1 
%  remodulate  and  send  it  to  Receiver  2  as  reference  signal. 

%  repeats  at  the  rest  of  the  receivers 
Rx  proced  =  zeros ( 4 , length (Rx  received));  %Buffer 
%  Begin  with  SIC,  then  perform  demodulation 
for  rev  count  =1:4  %  4 :  FOUR  receivers 

%  Interference  cancellation 

if  rev  count  ==  1  %  At  Receiver  One,  no  SIC  and  detection 
Rx  proced (1,:)=  Rx  received ( 1 ,:) ;  %  save  ref  signal  1 
elseif  rev  count  ==  2  %  Perform  RSIC  at  Receiver  Two 

Rx_proced (2 , : ) =Rx_received (2 , : ) - (beta!2*Rx_proced ( 1 , : ) ) ; 
elseif  rev  count  ==  3  %  At  Receiver  Three 
Rx_proced (3, : ) =Rx_received (3, : ) - (beta23*Rx_proced (2, : ) ) . . . 
-  (betal3  *  Rx_proced ( 1 , : ) ) ; 
elseif  rev  count  ==  4  %  At  Receiver  Four 
Rx_proced ( 4 ,  : ) =Rx_received ( 4 ,  : ) - (beta34  *Rx_proced ( 3 , : ) )  ... 
- (beta24*Rx_proced (2, : ) ) - (betal4*Rx_proced (1, : ) ) ; 

end 

if  ch  ==  2  %  OFDM  simulation 
%  OFDM  Demodulation 
for  n  =  1  :  n  symbols; 

r  ofdm_cyc  =  Rx  proced (:, (n-1) * (n  fft+n  cyc)+l  :... 

n* (n  fft+n  eye) ) ;  %  retrieve  each  ofdm  symbol 

r  ofdm  =  r  ofdm  eye ( : ,  n  cyc+l:n  fft+n  eye)  ';  %Remove  CP 
r  qpsk=fft(r  ofdm,n  fft)/sqrt(n  fft)  ;  %Normalized  FFT 
QPSK  ref ( : ,  (n-l)*n  fft+l:n*n  fft)  =  r  qpsk'; 
end 

Rx  proced  =  QPSK  ref; 

end 
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%  QPSK  detection  on  received  signal 

%  Need  to  flip  I  and  Q  here  to  have  correct  sequence! 

Q (rcv_count, : )  =  (real (Rx_proced (rcv_count, : ) ) <0) ; 

I (rcv_count, : )  =  (imag (Rx_proced (rcv_count,  : ) ) <0) ; 

%  Remodulate  for  transmition  (as  referecnce  signal) 
Rx_proced (rcv_count,  :)  =  ... 

( (I (rcv_count, : ) ==0) . * (Q (rcv_count, : ) ==0) * (exp (li*pi/ 4) ) . . . 

+ (I (rcv_count, : ) ==0) . * (Q (rcv_count, : ) ==1) * (exp (3*li*pi/ 4) ) ... 

+ (I (rcv_count, : ) ==1) .  *  (Q (rcv_count, : ) ==1) * (exp (5*li*pi/4) ) . . . 

+  (I (rcv_count,  : ) ==1)  . * (Q (rcv_count,  : ) ==0) * (exp  (7*li*pi/ 4) ) ) ; 

if  ch  ==2  %  OFDM  simulation 
%  OFDM  re-Modulation 
for  n  =  1  :  n  symbols; 

s  ofdm=sqrt (n  fft)*ifft (Rx  proced ( : , (n-1) *n  f ft+1 : . . . 
n  *  n  fft)',n  fft); 

s_ofdm_cyc  =  sqrt (n_fft/ (n_f ft+n_cyc) ) * . . . 

[s  ofdm(n  fft-n  cyc+l:n  fft, :);s  ofdm]; 

OFDM  ref ( : , (n-1 ) * (n  fft+n  eye) +1 : n* (n  fft+n  cyc))=... 
s_ofdm_cyc ' ; 

end 

Rx  proced  =  OFDM  ref; 
end  %  End  of  OFDM  simulation 
end  %  End  of  Demodulated-RSIC 

end  %  End  of  three  RSIC  techniques  loop 

0,0,  _ 

o  o 

%%  Final  demodulation  at  receiver  for  BER  calculation 
if  ch  ==  2  %  OFDM  simulation 
%  OFDM  Demodulation 
for  n  =  1  :  n  symbols; 

r  ofdm  cyc=Rx  proced (:, (n-1 )* (n  fft+n  cyc)+l:n*(n  fft+n  eye)); 
r  ofdm  =  r  ofdm  eye (  :  ,  n  cyc+l:n  fft+n  eye)  ';  %Remove  CP 
r  qpsk=fft(r  ofdm,n  fft) /sqrt (n  fft);  %Normalized  FFT 
QPSK  Rev ( : , (n-l)*n  fft+l:n*n  fft)  =  r  qpsk'; 
end 

Rx  proced  =  QPSK  Rev; 
end 

%  QPSK  Demodulation 

I_rx  =  ( real (Rx_proced) <0 ) ;  %Detection 
Q  rx  =  (imag(Rx  proced) <0); 

uncoded  bits  rx  =  zeros (4, 2*length (Rx  proced)); 
uncoded  bits  rx ( :  ,  1 : 2 : end)  =  Q  rx;  %converge  split  streams 

uncoded  bits  rx ( : , 2 : 2 : end)  =  I  rx;  ^converge  split  streams 

%%  Sum  up  bit  errors 

diff  =  uncoded  bits  -  uncoded  bits  rx;  %  find  the  errors 
Tx  Errors  =  Tx  Errors  +  sum (abs (dif f ) , 2 ) ;  %  2  :  sum  columns 
Tx  bits  =  Tx  bits  +  length (uncoded  bits); 
end  %  end  of  while  loop (keep  running  until  #  of  errors) 

%%  Calculate  Bit  Error  Rate  of  each  SNR 

BER(:,snr  count)  =  Tx  Errors. /Tx  bits;  %Calculate  BER  for  each  SNR 
end  %  end  of  each  SNR  calculation  loop 

o, _ 

o 

%%  save  BERs  of  each  SNR 
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%  save  BER  (all  four  channels)  of  ND-RSIC  and  D-RSIC  separately, 
if  RSIC_count  ==  1  %  Non-RSIC 

BER_Non_RISC  =  BER;  %  Save  BER  of  all  4  signals  (Non-RSIC) 
elseif  RSIC_count  ==  2  %  Non-Demodulated  RSIC 

BER_ND_RSIC  =  BER;  %save  BER  of  all  4  signals  (ND-RSIC) 
elseif  RSIC_count  ==  3  %  Demodulated  RSIC 

BER_D_RSIC  =  BER;  %  save  BER  of  all  4  signals  (D-RSIC) 
end  %  end  of  saving  BER  of  all  four  signals  in  three  RSIC 

o. _ 

o 

end  %  end  of  three  RSIC  techniques  loop 
%%  save  Receiver  2  BER  of  each  alphal2 

BER_ND_RSIC_Rcv2_alpha (alpha_count, : )  =  BER_ND_RSIC (2, : ) ;  %ND_RSIC 
BER_D_RSIC_Rcv2_alpha (alpha_count, : )  =  BER_D_RSIC (2 , : ) ;  %D_RSIC 

%%  Comparison  of  three  RSIC,  Show  Receiver  2 (same  as  two  receivers) 
figure ( 1 ) 

subplot (2,2, alpha_count) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erfc (sqrt (10 . A (EbNo/10) ) ) ; 
semi logy ( SNR, theoryBer,  *  — k ' ,  ' LineWidth ' , 2 ) ; 
hold  on 

%  Simulation  results 

semi logy ( SNR, BER  Non  RISC (2  ,  : )  ,  ' m-o ' ,  ' LineWidth '  ,  2 ) 
semilogy (SNR, BER~ND_RSIC (2,  : ) ,  'b-d'  ,  'LineWidth' ,2) 
semi logy (SNR, BER_D_RSIC (2,  : ) ,  'r-A '  ,  'LineWidth' ,2) 
xlabel ( ' SNR  (dB) ' ) ; 
ylabel ( ' BER ' ) ; 

legend( 'QPSK  Theory' ,  'Non-RSIC',  'ND-RSIC',  'D-RSIC'); 
if  ch  ==  1  %  QPSK  simulation 

title ([' Three  RSICs  in  a  two  QPSK  receiver  system',... 
sprintf ( ' \n ' )  ,  ' \alpha  { 12 }  =  ' ,  num2str (alphal2 (alpha  count) ) ,  .  .  . 
';  \gamma  {12}  =  ',  num2str (gammal2 ) ] ) ; 
elseif  ch  ==  2  %  OFDM  simulation 

title ([' Three  RSICs  in  a  two  OFDM  receiver  system',... 

sprintf (' \n '),  ' Alpha  { 12 }  =  ' ,  num2str (alphal2 (alpha  count) ) ,  .  .  . 

';  \gamma  {12}  =  ',  num2str (gammal2 ) ] ) ; 

end 

grid  on; 

%%  plot  ND-RSIC  all  four  signals  (in  one  given  alpha  value) 
figure  (2 ) 

subplot (2,2, alpha_count) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erfc (sqrt (10 . A (EbNo/10) )) ; 
semi logy ( SNR, theoryBer,  ' -k '  ,  ' LineWidth ' , 2 ) 
hold  on 

%  Simulation  results  (Non-Demudulated  RSIC) 
semilogy (SNR, BER_ND_RSIC (1, : ) , '-or' , 'LineWidth' ,2) 
semi logy (SNR, BER_ND_RSIC (2,  : ) ,  ' -db' ,  'LineWidth'  ,2) 
semilogy (SNR, BER_ND_RSIC (3,  : ) ,  '-Ag'  ,  'LineWidth' ,2) 
semilogy (SNR, BER_ND_RSIC (4, : ) , ' -*m' , 'LineWidth' ,2) 
xlabel ( ' SNR  (dB) ' ) ; 
ylabel ( ' BER ' ) ; 

legend ( 'QPSK  Theory ' ,  ' Rx  Signal  1 ' ,  ' Rx  Signal  2 ' ,  ' Rx  Signal  3 '  ,  .  .  . 
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'Rx  Signal  4  '  )  ; 
if  ch  ==  1  %  QPSK  simulation 

title (['Four  ND-RSIC  QPSK  receivers  demodulating  signals',... 
sprintf ( ' \n ' ) , ' \alpha_{ 12 }  =  '  ,  num2str (alphal2 (alpha_count) ) , . . . 
';  \gamma  {12}  =  ',  num2str (gammal2 ) ] ) ; 
elseif  ch  ==  2  %  OFDM  simulation 

title (['Four  ND-RSIC  OFDM  receivers  demodulating  signals',... 
sprintf ( ' \n ' )  ,  ' \alpha  { 12 }  =  ' ,  num2str (alphal2 (alpha  count) ) ,  .  .  . 
';  \ gamma  {12}  =  ',  num2str (gamma 12 )]) ; 

end 

grid  on 
hold  off 

%%  plot  Demodulated  RSIC  all  four  signals  (in  one  given  alpha  value) 
figure  ( 3 ) 

subplot (2,2, alpha_count) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erfc (sqrt (10 . A (EbNo/10) ) ) ; 
semi logy ( SNR, theoryBer, ' -k ' , ' LineWidth ' , 2 ) ; 
hold  on 

%  Simulation  results  (Demodulated  RSIC) 

semi logy ( SNR,  BER  D  RSIC ( 1 ,  : )  ,  ' -or ' ,  ' LineWidth '  ,  2 ) 

semi logy ( SNR, BER  D  RSIC (2  ,  : )  ,  ' -db ' ,  ' LineWidth ' , 2 ) 

semi logy ( SNR, BER  D  RSIC(3,:),'-Ag',' LineWidth ' , 2 ) 

semi logy (SNR, BER~D~RSIC (4, : ) , ' -*m' , 'LineWidth' ,2) 

xlabel ( ' SNR  (dB) ' ) ; 

ylabel ( ' BER ' ) ; 

legend ( 'QPSK  Theory ' ,  ' Rx  Signal  1 ' ,  ' Rx  Signal  2 ' ,  ' Rx  Signal  3  '  ,  .  .  . 

' Rx  Signal  4 ' )  ; 

%axis (  [  8  13  10A  (-6)  1.0]) 
if  ch  ==  1  %  QPSK  simulation 

title (['Four  ND-RSIC  QPSK  receivers  demodulating  signals',... 
sprintf ( ' \n ' ) ,  ' \alpha_{ 12 }  =  ' ,  num2str (alphal2 (alpha_count) )  ,  .  .  . 
';  \gamma  {12}  =  ',  num2str (gammal2 ) ] ) ; 
elseif  ch  ==  2  %  OFDM  simulation 

title (['Four  ND-RSIC  OFDM  receivers  demodulating  signals',... 
sprintf ( ' \n ' )  ,  ' \alpha  { 12 }  =  ' ,  num2str (alphal2 (alpha  count) ) ,  .  .  . 
';  \ gamma  {12}  =  ',  num2str (gamma 12 )]) ; 

end 

grid  on 
hold  off 

end  %  end  of  various  Alphal2  loop 

%%  Plot  BER  Vs.  SNR(dB)  Curve  on  logarithmic  scale 
%%  Non-Demodulated  RSIC  Receiver  2  in  various  alpha 
%%  Non-Demodulated  RSIC  Receiver  2  in  various  alpha 
figure ( 4 ) 
subplot  (1,2, 1 ) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erf c ( sqrt ( 1 0 . A (EbNo/ 1 0 ) ) ) ; 
semi logy ( SNR, theoryBer,  ' -k  '  ,  ' LineWidth ' , 2 ) ; 
hold  on 

%  Simulation  results  (Non-Demudulated  RSIC) 

semilogy ( SNR, BER  ND  RSIC_Rcv2  alpha ( 1 ,:),' -or ',' LineWidth ', 2 ) 
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semilogy ( SNR, BER  ND  RSIC_Rcv2  alpha (2 , : ) , ' -db ' , ' LineWidth ' , 2 ) 
semilogy ( SNR, BER  ND  RSIC  Rcv2  alpha ( 3 -Ag LineWidth ', 2 ) 
semilogy ( SNR, BER  ND  RSIC  Rcv2  alpha ( 4 -*m LineWidth ',  2 ) 
xlabel ( ' SNR  (dB)T);_ 
ylabel ( ' BER ' ) ; 

if  ch  ==  1  %  QPSK  simulation 

title ([' ND-RSIC  QPSK  Receiver  2  demodulation  with  various  alpha',, 
sprintf ( ' \n ' ) ,  '  (  \alpha  { 12 }  =  0.5  to  2  and  \gamma  { 12 }  =  '  ,  .  . 
num2str (gamma 12 ),')']); 
elseif  ch  ==  2  %  OFDM  simulation 

title ([' ND-RSIC  OFDM  Receiver  2  demodulation  with  various 
alpha ' , . . . 

sprintf ( ' \n ' )  ,  '  (  \alpha  { 12 }  =  0.5  to  2  and  \gamma  { 12 }  =  '  ,  .  . 
num2str (gammal2 ),')']); 

end 

legend('QPSK  Theory',  '\alpha_{12}  =  0.5',  '\alpha_{12}  =  1.0',... 

' \alpha_{ 12 }  =  1.5',  '\alpha_{12}  =2') 

axis ( [6  13  10A  (-6)  1.0] ) 
grid  on; 

%%  Demodulated  RSIC  Receiver  2  in  various  alpha 
figure  ( 4 ) 
subplot  (1,2,2) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erf c ( sqrt ( 1 0 . A (EbNo/ 1 0 ) ) ) ; 
semilogy ( SNR, theoryBer,  ' -k  '  ,  ' LineWidth ' , 2 ) ; 
hold  on 

%  Simulation  results  (Demudulated  RSIC) 

semilogy ( SNR,  BER  D  RSIC  Rcv2  alpha ( 1 ,:),' -or ',' LineWidth ', 2 ) 
semilogy ( SNR, BER  D  RSIC  Rcv2  alpha (2 ,:),' -db ',' LineWidth ', 2 ) 
semilogy ( SNR, BER  D  RSIC  Rcv2  alpha ( 3 ,:),' -Ag ',' LineWidth ', 2 ) 
semilogy ( SNR, BER  D  RSIC  Rcv2  alpha ( 4 ,:),' -*m ',' LineWidth ', 2 ) 
xlabel ('SNR  (dB)T)7 
ylabel ( ' BER ' ) ; 

if  ch  ==  1  %  QPSK  simulation 

title ([' D-RSIC  QPSK  Receiver  2  demodulation  with  various  alpha',., 
sprintf ( ' \n ' ) ,  '  (  \alpha  { 12 }  =  0.5  to  2  and  \gamma  { 12 }  =  '  ,  .  . 
num2str (gammal2 ),')']); 
elseif  ch  ==  2  %  OFDM  simulation 

title ([' D-RSIC  OFDM  Receiver  2  demodulation  with  various  alpha 
sprintf ( ' \n ' ) ,  '  (  \alpha  { 12 }  =  0.5  to  2  and  \gamma  { 12 }  =  '  ,  . . 
num2str (gammal2 ) , ' ) ' ] ) ; 

end 

legend('QPSK  Theory',  '\alpha_{12}  =  0.5',  '\alpha_{12}  =  1.0',... 

' \alpha_{ 12 }  =  1.5',  ’\alpha_{12}  =2') 

axis ([6  13  10A  (-6)  10A  (0) ] ) 
grid  on; 
toe 

%  ********************END  of  Sfmulat;j_on************************ 


43 


Figures  9  to  13  (QPSK  simulation  with  various  y  )  and  Figures  18  to  22  (OFDM 


simulation  with  various  y  )  are  generated  by  the  following  Matlab  codes: 


Q~'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k~k'k'k'k'k'k'k 

o 

%  Description: 

%  1.  Four  QPSK/OFDM  emitters  and  receivers 

%  2.  Three  techniques:  non-RSIC,  Non-Demodulated  and  Demodulated  RSIC 
%  3.  At  Receiver  2,  alphal2  varies  and  Gammal2  =  -0.05 (fixed) 

o 


clc;  clear  all;  close  all 
format  long; 

%  choose  simulation  model:  (l)QPSK;  (2) OFDM 

disp('RSIC  simulation  with  fixed  alpha  and  various  gamma  value') 
ch  =  input (' Choose  (l)QPSK  Simulation  or  (2) OFDM  Simulation  :  '); 

ch  gamma  =  input (' Choose  (1) Positive  or  (2) Negative  Gamma  Value:  '); 
tic  %  start  clock  to  measure  performance 
%  Bit  length 


if  ch  ==  1 


QPSK  simulaiton 


bit  length  =  le6; 
elseif  ch  ==  2  %  OFDM  simulation 
%  Bit  Length 
M=4 ; 

k=log2 (M) ; 
n  fft=64; 
n  cyc=16; 


%  Number  of  symbols  in  the  constellation 
i  Bits  per  symbol 

%  Size  of  FFT (number  of  OFDM  carriers) 

%  Length  of  cyclic  prefix  (OFDM) 


n  symbols  =  le4;  %  Number  of  symbols;n  symbols=(bit  count/k) /n  fft 
bit  length=k*n  fft*n  symbols;  %bl= (bits/symbol) * (FFTsize) * (symbols) 


end 

%  Range  of  SNR(dB)  to  simulate 
EbNo  =  -3:  1:  10; 

%  Convert  Eb/No  values  to  channel  SNR 
SNR  =  EbNo  +  10*logl0 (2 ) ; 


%%  Define 

variables  - 

interference 

amplitude 

gain 

%  Alpha 

alphal2  = 

.5;  % 

Signal 

1 

leaks 

to  Receiver  2 

alphal3  = 

.9;  % 

Signal 

1 

leaks 

into 

Receiver 

3 

alpha23  = 

.9;  % 

Signal 

2 

leaks 

into 

Receiver 

3 

alphal4  = 

.8;  % 

Signal 

1 

leaks 

into 

Receiver 

4 

alpha24  = 

.8;  % 

Signal 

2 

leaks 

into 

Receiver 

4 

alpha34  = 

.7;  % 

Signal 

3 

leaks 

into 

Receiver 

4 

%  Beta:  Reference  signal  amplitude  received  at  subsequent  receiver 


betal3  = 

LO 

Oh 

%Receiver 

1 

pass 

reference 

signal 

to 

receiver 

3 

betal4  = 

.95; 

^Receiver 

1 

pass 

reference 

signal 

to 

receiver 

4 

beta23  = 

.9; 

%Receiver 

2 

pass 

reference 

signal 

to 

receiver 

3 

beta24  = 

.9; 

%Receiver 

2 

pass 

reference 

signal 

to 

receiver 

4 

beta34  = 

.  65; 

%Receiver 

3 

pass 

reference 

signal 

to 

receiver 

4 

%  Various  gammal2  (Gamma  =  Alpha  -  Beta) 
if  ch  gamma  ==  1  %  positive  gamma  value 

gammal2  =  [0.01  0.05  0.1  0.15] ; 
elseif  ch  gamma  ==  2  %  negative  gamma  value 
gammal2  =  -0.1:  0.025:  -0.025; 
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end 

%%  create  BER  buffers  (to  increase  Matlab  speed,  not  necessary) 

BER  =  zeros ( 4 , length ( SNR) ) ;  %  4:  FOUR  receivers 

%save  Rcvr2  SNR  of  various  gammal2 

BER  ND  RSIC  Rcv2  gamma=  zeros ( length (gammal2 ), length ( SNR) ) ; 

BER  D  RSIC_Rcv2  gamma  =  zeros ( length (gammal2 ), length ( SNR) ) ; 

%%  Start  various  Gammal2  calculation  loop 

for  gamma  count  =1:  1:  length (gammal2 )  %  Calculate  each  gammal2 

betal2  =  alphal2  -  gammal2 (gamma  count) ;  %  various  gammal2 
for  RSIC  count  =1  :  1  : 3  %  Do  three  RSIC  techniques  one  at  a  time 
%  1:  Non-RSIC;  2:  Non-demodulated  RSIC;  3:  Demodulated  RSIC 
%%  Start  the  main  calculation  loop  for  each  SNR 
for  snr  count  =1:  1:  length (SNR) 

%  Initiate  values  for  error  buffer 

Tx  Errors  =  [0;  0;  0;  0];  %  FOUR  rows  =  FOUR  receivers 
Tx_bits  =  [0;  0;  0;  0] ;  %  FOUR  rows  =  FOUR  receivers 

%  Keep  running  until  we  get  100  errors  in  Signal  1 
while  Tx  Errors (1,:)  <  100 

%%  Generate  information  bits;  4  rows  for  4  signals (emitters) 
uncoded_Jbits  =  round  (rand  ( 4  ,  bit_length)  )  ; 

%%  Split  each  signal  into  two  streams,  for  Quadrature  Carriers 
I  =  uncoded__bits  ( : ,  1 :  2  :  end)  ;  %  get  odd  bits 

Q  =  uncoded_bits ( : , 2 ; 2 : end) ;  %  get  even  bits 

%  QPSK  modulator  set  to  pi/4  radians  constellation (emitter) 

Tx  sig  = . . . 

(  (I==0)  . * (Q==0) * (exp (li*pi/4) )  . . . 

+  (I==0)  . * (Q==l ) * (exp (3*li*pi/4) )  .  .  . 

+  (I==l)  . * (Q==l) * (exp (5*li*pi/4) )  .  .  . 

+  (1=1)  .  *  (Q==0)  *  (exp  (7*li*pi/4)  )  )  ; 
if  ch  ==  2  %  OFDM  simulation 
%%  OFDM  Modulation 
for  n  =  1  :  n  symbols; 

s  of dm  =  sqrt (n  f f t) *if f t (Tx  sig ( : , (n-1) *n  f ft+1 : . . . 

n*n  fft) ' ,n  fft) ;  %Normalized  IFFT  operation 
s  ofdm  eye  =  sqrt (n  fft/ (n  fft+n  eye))  *  [s  ofdm(n  fft... 

-n  cyc+1  :  n  fft, :);  s  ofdm];  %Appending  cyclic  prefix 
OFDM  tx ( : , (n-1) * (n  fft+n  eye) +1  :  n* (n  fft+n  eye) ) . . . 

=  s  ofdm  eye';  %  Combine  subcarriers  into  one  to  tx 

end 

Tx_sig  =  OFDM_tx; 

end 


%%  Gaussian  Link  Noise 

No  =  1/10 A ( SNR ( snr  count) /10);  %  Noise  variance 

Noise  =  sqrt (No/2 )*.. .  %  4  ramdom  noises  for  four  channels 

(randn (4, length (Tx  sig) ) +li*randn (4, length (Tx  sig))); 

%%  Receiver  intercepted  signals 

%  including  target  SOI  and  interferences  from  adjacent  emitters 
Rx  received  =  zeros (4, length (Tx  sig));  %  Create  buffer 
%  Receiver  1:  Signal  1  +  Noise  1 
Rx_received ( 1 , : )  =  Tx_sig(l,:)+  Noise ( 1, :); 

%  Receiver  2:  Signal  2  +  (alphal2  *  SI)  +  Noise  2 
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Rx  received (2 ,: )  =  Tx  sig (2, :)+... 

(alphal2  *  Tx_sig ( 1 , : ) ) +  Noise(2,:); 

%  Receiver  3:  Signal  3  +(alpha23  *  S2)+(alphal3  *  Sl)+Noise  3 
Rx_received ( 3 , : )  =  Tx_sig(3,:)+  (alpha23  *  Tx_sig (2 . 

(alphal3  *  Tx_sig (1 , : ) ) +  Noise(3,:); 

%  Receiver  4:  S4  +(a34  *  S3)  +(a24  *  S2)  +(al4  *  SI)  +Noise  4 
Rx  received ( 4 ,: )  =  Tx  sig(4,:)+  (alpha34  *  Tx  sig ( 3 ,:))  +  .. . 

(alpha24  *  Tx  sig (2, : ) ) +  (alpha 14  *  Tx  sig (1, : ) ) +  Noise  (4,  : ) ; 


%%  Reference-based  successive  interference  cancellation 
%  run  one  RSIC  technique  at  a  time 
%%  Non-RSIC 
if  RSIC_count  ==  1 

Rx  proced  =  Rx  received;  %  no  cancellation  applied 

%%  Non-Demodulated  RSIC 
elseif  RSIC  count  ==  2 

%  Signal  1  is  used  as  reference  directly;  Intercepted  signals 
%  will  be  subtract  by  reference  signal  with  a  factor  Beta 
Rx  proced  =  Rx  received;  %  Create  buffer  and  save  Signal  1 
Rx_proced (2 , : ) =Rx_received (2 , : ) - (betal2  *  Rx_proced ( 1 , : ) ) ;  %S2 
Rx_proced (3, : ) =Rx_received (3, : ) - (beta23  *  Rx_proced (2 ,:)).. . 

-  (betal3  *  Rx_proced ( 1 , : ) ) ;  %  get  Signal  3  only 

Rx_proced ( 4 , : )  =  Rx_received ( 4 , : )  -  (beta34 *Rx_proced ( 3 ,:)).. . 

-  (beta24  *  Rx_proced (2 , : ) )  -  (betal4*Rx_proced (1, : ) ) ;  %  S4 

%%  Demodulated  RSIC 
elseif  RSIC  count  ==  3 

%  Demodulate  1st  received  signal  at  Receiver  One, 

%  remodulate  and  send  it  to  Receiver  Two  as  reference  signal. 

%  repeats  at  the  rest  of  the  receivers 

Rx  proced  =  zeros ( 4 , length (Rx  received));  %Buffer;  4  signals 

%  Begin  with  SIC,  then  perform  demodulation 
for  rev  count  =1:4  %  4 :  FOUR  receivers 

%  Interference  cancellation 

if  rev  count  ==  1  %  At  Receiver  One,  no  SIC  and  detection 
Rx  proced (1,:)=  Rx  received ( 1 ,:) ;  %  save  RSl(tNoisel) 
elseif  rev  count  ==  2  %  Perform  RSIC  at  Receiver  Two 

Rx_proced (2, : ) =Rx_received (2, : ) - (betal2 *Rx_proced (1, : ) ) ; 
elseif  rev  count  ==  3  %  At  Receiver  Three 

Rx_proced (3, : ) =  Rx_received (3, : ) - (beta23*Rx_proced (2, : ) ) . . . 

-  (betal3  *  Rx_proced ( 1 , : ) ) ; 
elseif  rev  count  ==  4  %  At  Receiver  Four 
Rx_proced ( 4 ,  : ) =Rx_received ( 4 ,  : ) - (beta34 *Rx_proced ( 3 ,:))..  . 
- (beta24*Rx_proced (2, : ) ) - (beta!4*Rx_proced (1, : ) ) ; 

end 

if  ch  ==  2  %  OFDM  simulation 
%  OFDM  Demodulation 
for  n  =  1  :  n  symbols; 

r  ofdm  eye  =  Rx  proced ( : , (n-1) * (n  fft+n  eye) +1  : . . . 

n* (n  fft+n  eye) ) ;  %  retrieve  each  ofdm  symbol 

r  ofdm  =  r  ofdm  eye ( : , n  cyc+l:n  fft+n  eye) ';  %Remove  CP 
r  qpsk=fft(r  ofdm, n  fft)/sqrt(n  fft) ;  %Normalized  FFT 
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QPSK  ref ( :  ,  (n-l)*n  fft+l:n*n  fft)  =  r  qpsk'; 
end 

Rx  proced  =  QPSK  ref; 

end 

%  QPSK  detection  on  received  signal 

%  Need  to  flip  I  and  Q  here  to  have  correct  sequence! ! 

Q (rcv_count, : )  =  (real(Rx_proced(rcv_count,:))<0); 

I (rcv_count, : )  =  (imag (Rx_proced (rcv_count, : ) ) <0) ; 

%  Remodulate  for  transmition  (as  referecnce  signal) 
Rx_proced (rcv_count,  :)  =  ... 

( (I (rcv_count,  : ) ==0)  . * (Q (rcv_count,  : ) ==0) * (exp (li*pi/ 4) )  .  .  . 

+  (I  (rcv_count,  :  )  =0)  .  *  (Q  (rcv_count,  :  )  ==1)  *  (exp  (3*li*pi/4)  )  .  .  . 

+ (I (rcv_count, :)==!) . * (Q (rcv_count, : ) ==1) * (exp (5*li*pi/4) ) . . . 

+  (I (rcv_count,  : ) ==1)  . * (Q (rcv_count,  : ) ==0) * (exp  (7*li*pi/ 4) ) ) ; 

if  ch  ==2  %  OFDM  simulation 
%  OFDM  re-Modulation 
for  n  =  1  :  n  symbols; 

s  ofdm  =  sqrt (n  fft)*ifft(Rx  proced (: , (n- 

1) *n_fft+l: . . . 

n*n  f f t )  ' , n  fft); 

s_ofdm_cyc  =  sqrt (n_f ft/ (n_f ft+n_cyc) ) * . . . 

[s  ofdm(n  fft-n  cyc+l:n  fft,:);s  ofdm]; 

OFDM  ref ( : , (n-l)*(n  fft+n  cyc)+l:n*(n  fft  + 

n_cyc) ) = . . . 

s_ofdm_cyc ' ; 

end 

Rx  proced  =  OFDM  ref; 
end  %  End  of  OFDM  simulation 
end  %  End  of  Demodulated-RSIC 

end  %  End  of  three  RSIC  techniques  loop 

0,0,  _ 

o  o 

%%  Final  demodulation  at  receiver  for  BER  calculation 

if  ch  ==  2  %  OFDM  simulation 
%  OFDM  Demodulation 
for  n  =  1  :  n  symbols; 

r  ofdm  cyc=Rx  proced (:,  (n-1 )* (n  fft+n  cyc)+l:n*(n  fft+n  eye)) 
r  ofdm  =  r  ofdm  eye ( : , n_cyc+l : n  fft+n_cyc) ';  %Remove  CP 
r  qpsk=fft(r  ofdm,n  fft) /sqrt (n  fft);  %Normalized  FFT 
QPSK  Rev ( : , (n-l)*n  fft+l:n*n  fft)  =  r  qpsk'; 
end 

Rx  proced  =  QPSK  Rev; 

end 

%  QPSK  Demodulation 

I  rx  =  (real (Rx  proced) <0);  %Detection 
Q  rx  =  (imag(Rx  proced) <0); 

uncoded  bits  rx  =  zeros (4, 2*length (Rx  proced)); 

uncoded  bits  rx ( : ,  1 : 2 : end)  =  Q  rx;  ^converge  split  streams 
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uncoded  bits  rx ( : ,  2 : 2 : end)  =  I  rx;  %converge  split  streams 
%%  Sum  up  bit  errors 

diff  =  uncoded  bits  -  uncoded  bits  rx;  %  find  the  errors 
Tx  Errors  =  Tx  Errors  +  sum (abs (dif f ) , 2 ) ;  %  2  :  sum  columns 
Tx  bits  =  Tx  bits  +  length (uncoded  bits); 
end  %  end  of  while  loop 

%%  Calculate  Bit  Error  Rate  of  each  SNR 

BER(:,snr  count)  =  Tx  Errors. /Tx  bits;  %Calculate  BER  for  each  SNR 
end  %  end  of  each  SNR  calculation  loop 

o, _ 

o 

%%  save  BERs  of  each  SNR 

%  save  BER  (all  four  channels)  of  ND-RSIC  and  D-RSIC  separately, 
if  RSIC_count  ==  1  %  Non-RSIC 

BER_Non_RISC  =  BER;  %  Save  BER  of  all  4  signals  (Non-RSIC) 
elseif  RSIC  count  ==  2  %  Non-Demodulated  RSIC 

BER_ND_RSIC  =  BER; %  save  BER  of  all  4  signals (ND-RSIC) 
elseif  RSIC  count  ==  3  %  Demodulated  RSIC 

BER_D_RSIC  =  BER;  %  save  BER  of  all  4  signals  (D-RSIC) 
end  %  end  of  saving  BER  of  all  four  signals  in  three  RSIC. 

o, _ 

o 

end  %  end  of  three  RSIC  techniques  loop 
%%  save  Receiver  2  BER  of  each  alphal2 

BER_ND_RSIC_Rcv2_gamma (gamma_count, : )  =  BER_ND_RSIC (2 , : ) ;  %ND_RSIC 
BER_D_RSIC_Rcv2_gamma (gamma_count, : )  =  BER_D_RSIC (2 , : ) ;  %D_RSIC 

%  Plot  Figures 

%%  Comparison  of  three  RSIC,  show  Receiver  2 
figure  ( 1 ) 

subplot (2,2, gamma_count) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erfc (sqrt (10 . A (EbNo/10) ) ) ; 
semi logy ( SNR, theoryBer,  ' -k  '  ,  ' LineWidth ' , 2 ) ; 
hold  on 

%  Simulation  results 

semi logy ( SNR, BER  Non  RISC (2 , ; ) , 'm-o ' , ' LineWidth ' , 2 ) 
semi logy (SNR, BER~ND_RSIC (2,  : ) ,  'b-d' ,  'LineWidth'  ,2) 
semi logy (SNR, BER_D_RSIC (2,  : ) ,  ' r-A '  ,  'LineWidth' ,2) 
xlabel ( ' SNR  (dB) ' ) ; 
ylabel ( ' BER ' ) ; 

legend( 'QPSK  Theory' ,  'Non-RSIC',  'ND-RSIC',  'D-RSIC'); 
if  ch  ==  1  %  QPSK  simulation 

title ([' Three  RSICs  in  a  two  QPSK  receiver  system',... 
sprint f ( ' \n ' )  ,  ' \ alpha  { 12 }  =  ' ,  num2str (alpha 12 ) ,  .  .  . 

';  \gamma  {12}  =  ',  num2str (gamma 12 (gamma  count))]); 
elseif  ch  ==  2  %  OFDM  simulation 

title ([' Three  RSICs  in  a  two  OFDM  receiver  system',... 
sprint f ( ' \n ' ) , ' Alpha  { 12 }  =  ' ,  num2str (alphal2 ) , . . . 

';  \gamma  {12}  =  ',  num2str (gammal2 (gamma  count))]); 

end 

grid  on; 

%%  plot  Non-Demodulated  RSIC  all  four  signals  (in  one  given  alpha 
value) 
figure (2 ) ; 
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subplot (2,2, gamma_count) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erfc (sqrt (10 . A (EbNo/10) ) ) ; 
semi logy ( SNR, theoryBer,  ' -k '  ,  ' LineWidth ' , 2 ) 
hold  on 

%  Simulation  results  (Non-Demudulated  RSIC) 
semi  logy  (SNR,  BER_ND_RSIC  ( 1 ,  :  )  ,  '  -or'  ,  'LineWidth'  ,2) 
semi logy (SNR, BER_ND_RSIC (2, : ) , ' -db' , 'LineWidth' ,2) 
semi logy (SNR, BER_ND_RSIC (3, : ) , ' -Ag' , 'LineWidth' ,2) 
semi logy (SNR, BER_ND_RSIC (4, : ) , '-*m' , 'LineWidth' ,2) 
xlabel ( ' SNR  (dB) ' ) ; 
ylabel ( ' BER ' ) ; 

legend ( ' QPSK  Theory ' ,  ' Rx  Signal  1 ' ,  ' Rx  Signal  2 ' ,  ' Rx  Signal  3  '  ,  .  .  . 

'Rx  Signal  4 ' )  ; 
if  ch  ==  1  %  QPSK  simulation 

title (['Four  ND-RSIC  QPSK  receivers  demodulating  signals',... 
sprint f ( ' \n ' ) , ' \ alpha  { 12 }  =  ' ,  num2str (alpha 12 ) , . . . 

';  \gamma  {12}  =  ',  num2str (gammal2 (gamma  count))]); 
elseif  ch  ==  2  %  OFDM  simulation 

title (['Four  ND-RSIC  OFDM  receivers  demodulating  signals',... 
sprintf ( ' \n ' ) , ' \alpha_{ 12 }  =  ',  num2str (alphal2 ) , . . . 

';  \gamma  {12}  =  ',  num2str (gammal2 (gamma  count))]); 
end 

grid  on 
hold  off 

%%  plot  Demudulated  RSIC  all  four  signals  (in  one  given  alpha  value) 
figure  ( 3 )  ; 

subplot (2 , 2 , gamma  count) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erfc (sqrt (10 . A (EbNo/10) )) ; 
semi logy ( SNR, theoryBer, ' -k ' , ' LineWidth ' , 2 ) ; 
hold  on 

%  Simulation  results  (Demudulated  RSIC) 

semi logy ( SNR,  BER  D  RSIC (1 ,  : )  ,  ' -or ' ,  ' LineWidth '  ,  2 ) 

semi logy ( SNR, BER  D  RSIC (2  ,  : )  ,  ' -db ' ,  ' LineWidth ' , 2 ) 

semi logy ( SNR, BER  D  RSIC(3,:),'-Ag',' LineWidth ' , 2 ) 

semi logy (SNR, BER_D_RSIC ( 4 , : ) , ' -*m' , 'LineWidth' ,2) 

xlabel ( ' SNR  (dB) ' ) ; 

ylabel ( ' BER ' ) ; 

legend ( 'QPSK  Theory ' ,  ' Rx  Signal  1 ' ,  ' Rx  Signal  2 ' ,  ' Rx  Signal  3  '  ,  .  .  . 

' Rx  Signal  4 ' )  ; 
if  ch  ==  1  %  QPSK  simulation 

title (['Four  D-RSIC  QPSK  receivers  demodulating  signals',... 
sprintf (' \n '),' \alpha_{ 12 }  =  ',  num2str (alphal2 ) , . . . 

';  \gamma  {12}  =  ',  num2str (gammal2 (gamma  count))]); 
elseif  ch  ==  2  %  OFDM  simulation 

title ([' Four  D-RSIC  OFDM  receivers  demodulating  signals',... 
sprintf (' \n '), ' \ alpha  { 12 }  =  ' ,  num2str (alpha 12 ) , . . . 

';  \gamma  {12}  =  ',  num2str (gamma 12 (gamma  count))]); 
end 

grid  on 
hold  off 

end  %  end  of  various  Gammal2  loop 


49 


-or '  , 

' LineWidth ' 

,2) 

-db '  , 

' LineWidth ' 

,2) 

-  Ag '  , 

' LineWidth ' 

,2) 

-*m'  , 

' LineWidth ' 

,2) 

%%  Plot  BER  Vs.  SNR(dB)  Curve  on  logarithmic  scale 
%%  Non-Demodulated  RSIC  Receiver  2  with  various  gamma 
figure ( 4 ) 
subplot  (1,2, 1 ) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erf c ( sqrt ( 1 0 . A (EbNo/ 1 0 ) ) ) ; 
semi logy ( SNR, theoryBer,  ' -k  '  ,  ' LineWidth '  ,  2  )  ; 
hold  on 

%  Simulation  results  (Non-Demudulated  RSIC) 
semilogy ( SNR, BER  ND  RSIC  Rcv2  gamma (1,: 
semilogy ( SNR, BER  ND  RSIC  Rcv2  gamma(2,: 
semilogy ( SNR, BER  ND  RSIC  Rcv2  gamma(3,: 
semilogy ( SNR, BER  ND  RSIC  Rcv2  gamma(4,: 
xlabel ( ' SNR  (dB)  '  )  ; 
ylabel ( ' BER ' ) ; 

if  ch  ==  1  %  QPSK  simulation 

title ([' ND-RSIC  QPSK  Receiver  2  BER  with  various  \gamma  {12} 
'  (  \alpha_{  12  }  =  '  ,  num2str  (alphal2  ),')']); 
elseif  ch  ==  2  %  OFDM  simulation 

title ([' ND-RSIC  OFDM  Receiver  2  BER  with  various  \gamma  {12} 
'  (  \alpha_{  12  }  =  '  ,  num2str  (alphal2  ),')']); 

end 

if  ch  gamma  ==  1  %  positive  gamma  value 

legend ('QPSK  Theory \gamma  {12}  =  0.01', 

' \ gamma  {12}  =  0.1',  ' \ gamma  {12}  = 

elseif  ch  gamma  ==  2  %  negative  gamma  value 

legend ('QPSK  Theory ',' \gamma  {12}  =  -0.1', 

0.075' ,.. . 

' \ gamma  {12}  =  -0.05',  ' \ gamma  {12} 

end 

axis ([6  13  10A  (-6)  10A  (-1) ] ) 
grid  on; 


' \gamma 
0.15') 


{12}  =  0. 


' \gamma  { 12 }  =  - 
=  -0.025’  ) 


%%  Demodulated  RSIC  Receiver  2  in  various  alpha 
figure  ( 4 ) 
subplot (1,2,2) 

%  Theoretical  BER 

theoryBer  =  0 . 5*erfc (sqrt (10 . A (EbNo/10) ) ) ; 
semilogy ( SNR, theoryBer,  ' -k  '  ,  ' LineWidth ' , 2 ) ; 
hold  on 

%  Simulation  results  (Demudulated  RSIC) 

semilogy ( SNR, BER  D  RSIC  Rcv2  gamma ( 1 ,:),' -or ',' LineWidth ', 2 ) 
semilogy ( SNR, BER  D  RSIC  Rcv2  gamma (2 ,:),' -db ',' LineWidth ', 2 ) 
semilogy ( SNR, BER  D  RSIC  Rcv2  gamma ( 3 ,:),' -Ag ',' LineWidth ', 2 ) 
semilogy ( SNR, BER  D  RSIC  Rcv2  gamma ( 4 ,:),' -*m ',' LineWidth ', 2 ) 
xlabel ('SNR  (dB)T)7 
ylabel ( ' BER ' ) ; 

if  ch  ==  1  %  QPSK  simulation 

title ([' D-RSIC  QPSK  Receiver  2  BER  with  various  \gamma  {12} 
'  ( \alpha_{ 12 }  =  ' , num2str (alpha 12 )  ,')']); 
elseif  ch  ==  2  %  OFDM  simulation 

title ([' D-RSIC  OFDM  Receiver  2  BER  with  various  \gamma  {12} 
'  ( \alpha_{ 12 }  =  ' , num2str (alpha 12 )  ,')']); 

end 


05  ' 
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if  ch  gamma  ==  1  %  positive  gamma  value 

legend('QPSK  Theory \gamma  {12}  =  0.01', 

' \ gamma  {12}  =  0.1',  ' \ gamma  {12}  = 

elseif  ch  gamma  ==  2  %  negative  gamma  value 

legend('QPSK  Theory ',' \gamma  {12}  =  -0.1', 
0.075'  ,..  . 

' \ gamma  {12}  =  -0.05',  ' \ gamma  {12} 

end 

axis ([6  13  10A  (-6)  10A  (-1) ] ) 

grid  on; 

toe 


' \gamma_{ 12 } 
0.15') 

' \ gamma  { 12 } 
=  -0.025') 


0.05'  ,..  . 


o, 

o 


'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k'k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k-k'k'k'k'k'k-k'k-k'k-k'k-k'k'k'k'k'k-k'k'k'k-k'k'k'k-k'k 
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